Highly Frustrated Magnetic Clusters: The kagome on a sphere 
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We present a detailed study of the low-energy excitations of two existing finite-size realizations 
of the planar kagome Heisenberg antiferromagnet on the sphere, the cuboctahedron and the icosi- 
dodecahedron. After highlighting a number of special spectral features (such as the presence of 
low-lying singlets below the first triplet and the existence of localized magnons) we focus on two 
major issues. The first concerns the nature of the excitations above the plateau phase at 1/3 of the 
saturation magnetization Ms. Our exact diagonalizations for the s — 1/2 icosidodecahedron reveal 
that the low-lying plateau states are adiabatically connected to the degenerate collinear "up-up- 
down" ground states of the Ising point, at the same time being well isolated from higher excitations. 
A complementary physical picture emerges from the derivation of an effective quantum dimer model 
which reveals the central role of the topology and the intrinsic spin s. We also give a prediction 
for the low energy excitations and thermodynamic properties of the spin s — 5/2 icosidodecahe- 
dron Mo72Fe3o. In the second part we focus on the low-energy spectra of the s > 1/2 Heisenberg 
model in view of interpreting the broad inelastic neutron scattering response reported for Mo72Fe3o. 
To this end we demonstrate the simultaneous presence of several broadened low-energy "towers of 
states" or "rotational bands" which arise from the large discrete spatial degeneracy of the classical 
ground states, a generic feature of highly frustrated clusters. This semiclassical interpretation is 
further corroborated by their striking symmetry pattern which is shown, by an independent group 
theoretical analysis, to be a characteristic fingerprint of the classical coplanar ground states. 
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PACS numbers: 75.50.Xx,75.10.Jm,75.40.Mg 



I. INTRODUCTION 

The field of highly frustrated magnetism has received 
a growing theoretical and experimental interest in re- 
cent yearsi^. One of the central motifs in the planar 
kagome and similarly frustrated Heisenberg antifcrro- 
magnets (AFM's) which readily differentiates them from 
unfrustrated (e.g. collinear) ones, is the proliferation of 
an extensive family of low-energy singlets below the low- 
est triplet excitation'^''*. One interpretation for the origin 
of these singlets has emerged from Resonating Valence 
Bond (RVB) type of arguments^ for the s — 1/2 kagome 
AFM. For higher spins, purely classical considerations as- 
sert that the singlets stem from the splitting by quantum 
fluctuations of the extensively degenerate family of Neel 
ordered (3-sublattice) ground stated. Both interpreta- 
tions rest on the notion of a local degeneracy which stems 
from the frustrated corner-sharing topology of these lat- 
tices. In this regard, it appears that the proliferation of 
singlets is only one particular manifestation of this local 
degeneracy since similarly dense low-energy excitations 
are manifested in the whole magnetization range. 

On the other hand, some understanding for the ground 
state itself has been established. Exact Diagonalization 
(ED) results suggest that the ground state of the s — 1/2 
kagome AFM is a disordered spin liquid with a very 
small spin gap'^'^ (if any). For s > 1/2, semi-classical 
approaches predict that an extensive subset of coplanar 
states is first selected in 1/s while the -\/3 x -\/3 ordered 
state is stabilized in higher orders through the order-by- 



disorder mechanism^i^. In a magnetic field, the ground 
state may exhibit a number of interesting phases. These 
include the presence of an extensively degenerate family 
of localized magnons which result in macroscopic mag- 
netization jumps at the saturation field, as well as the 
stabilization of spin gaps and the associated fractional 
magnetization plateaux. For a first understanding of the 
nature of these plateaux a perturbative expansion around 
the degenerate Ising point was first employed by Cabra et 
al.^ for the kagome. This approach was recently extended 
by Bergman et ali^i^ to other frustrated systems, such 
as the pyrochlore AFM. Here, the anisotropy terms are 
treated perturbatively, and the emerging splitting of the 
degenerate Ising manifold is effectively cast into a Quan- 
tum Dimer Model (QDM) on the dual lattice. 

At the same time, it is well known that some pre- 
cursors of the excitation spectra of frustrated and un- 
frustrated AFM's are already embodied in the spectra 
of small system sizes (see for instance Ref. [l2l ). It has 
come therefore with no surprise that a number of phe- 
nomena that are manifest in kagome-like AFM's have 
also emerged in the research field of highly frustrated 
nanomagnet o^'^'^'^'^^1^^'^^1^^'^^ . These are realizations of 
zero-dimensional molecular-size magnets which consist of 
a finite number of strongly interacting transition metal 
ions, with the isotropic Heisenberg exchange being the 
dominant energy term. Thus, in addition to their great 
relevance in the context of nanomagnctism and the grow- 
ing interest for potential applications in quantum com- 
puting^*^, information storageSi and magnetic imaging—. 



molecular nanomagnets can also provide a suitable plat- 
form for addressing theoretical questions and testing 
ideas from the more general context of frustrated mag- 
netism. 

In this work, we focus on two magnetic molecule real- 
izations of the Heisenberg kagome AFM on the sphere. 
The first consists of 8 corner-sharing triangles and is re- 
alized in the Cui2Lag^'^ cluster with 12 Cu^+ s — 1/2 
ions occupying the vertices of a symmetric cuboctahe- 
dron (cf.Fig. [T]). The spin topology of this cluster is 
identical to the 12-site kagome wrapped on a torus (cf. 
Fig. [T6|). The second cluster is one of the largest frus- 
trated molecules synthesized to date, namely the giant 
Keplerate Mo72Fe3o system^!. This features an array 
of thirty s = 5/2 Fe'^"'" ions occupying the vertices of 
twenty corner-sharing triangles spanning an almost per- 
fect icosidodecahedron (cf.Fig. [T|). Interestingly, its quan- 
tum s — 1/2 analogue, M072V30, consisting of V^"*" ions 
has also been synthesized quite recentl y^^i^^ . We may 
note here that the cuboctahedron and the icosidodecahe- 
dron can be thought of as two existing positive curvature 
(with n — A and 5 respectively) counterparts of Elser 
and Zeng's21 generalization of the kagome structure on 
the hyperbolic plane where each hexagon is replaced by 
a polygon of n sides with n > 6. 

Among the above highly frustrated clusters, Mo72Fe3o 
has been the most investigated so far, both theoreti- 
cally and experimentally. The exchange interactions in 
Mo72Fe3o are quite small, J/ks — 1.57 K^''., and this has 
allowed for the experimental observation of a M = Ms/ 3 
plateau at H ~ 5.9 Tesla which has been explained 
classically by Schroder et al}'' . In addition, this clus- 
ter manifests a very broad Inelastic Neutron Scattering 
(INS) response as shown by Garlea et al^. On the 
other hand, M072V30 has a much stronger AFM exchange 
J/ks — 250 K— 12^, and thus is not well suited for the 
observation of the field-induced plateau. However, its 
low-energy excitation spectrum can still be investigated 
by INS experiments (which, to our knowledge, have not 
been performed so far). As to the s = 1/2 cuboctahedron 
Cui2LagSl, we are not aware of any magnetic measure- 
ments reported so far on this cluster. 



The main magnetic properties of the present clusters 
can be explained very well by the isotropic Heisenberg 
model with a single AFM exchange parameter J, i.e. 





FIG. 1: (Color online) Schematic representation of the cuboc- 
tahedron (left) and the icosidodecahedron (right). The first 
consists of 12 vertices, 24 edges, 6 square and 8 triangular 
faces, while the latter consists of 30 vertices, 60 edges, 12 
pentagons and 20 corner-sharing triangles. 
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where, as usual, (ij) denotes pairs of mutually interact- 
ing spins s at sites i and j. Other terms such as single- 
ion anisotropy (for s > 1/2) or Dzyaloshinsky-Moriya 
interactions must be present as well in the present clus- 
ters, but they are expected to be much smaller than the 
exchange interactions and thus they can be neglected. 
Here, as a simple theoretical tool to understand some of 
the properties of the Heisenberg model, it will be very ex- 
pedient to introduce some fictitious exchange anisotropy, 
i.e. extend Eq. ([T]) to its more general XXZ variant 
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where J^y, Jz denote the transverse and longitudinal ex- 
change parameters respectively. In what follows we de- 
note a = Jxy/Jz- 

The main results presented in this article are of direct 
relevance to the experimental findings in Mo72Fe3o men- 
tioned above and thus span two major themes. The first 
deals with the nature of the low-lying excitations above 
the M — Ms/Z plateau phase. For the s — 1/2 icosi- 
dodecahedron we show that all these excitations are adi- 
abatically connected to coUinear "up-up-down" (hence- 
forth "uud") Ising ground states (GS's), at the same time 
being well isolated from higher levels by a relatively large 
energy gap. We argue that this feature must be spe- 
cial to the topology of the icosidodecahedron and that 
it must survive for s = 5/2 as well. This prediction 
can be verified experimentally by a measurement of the 
low-temperature specific heat and the associated entropy 
content at the plateau phase of Mo72Fe3o. A comple- 
mentary physical picture will emerge by performing a 
high order perturbative expansion in a, in the spirit of 
Refs. lol llOlilll . and by deriving and solving to lowest or- 
der the corresponding effective QDM on the dual clusters. 
The dependence of the model parameters on a and s is 
also found and given explicitly. 

Our second theme concerns the origin of the broad 
INS response reported for Mo72Fe3o^^. Previous theories 
based on the excitations of the rotational band mode P^i^^ 
or on spin wave calculations'^'^''^^ predict a small number 
of discrete excitation lines at low temperatures and thus 
cannot explain the broad INS features. Our interpreta- 
tion of this behavior is based on the notion of the simul- 
taneous presence of several rotational bands or towers of 
states at low energies which originate from the large de- 
gree of classical degeneracy, a generic feature of highly 



frustrated systems. Indeed, our exact diagonalizations 
demonstrate the existence of an unusually high density 
of low-energy excitations manifesting in the full magne- 
tization range. A detailed group theory analysis reveals 
that the low-energy spectra are of semiclassical origin up 
to a relatively large energy cutoff. We will also show 
that the symmetry of the corresponding excitations for 
s — 1/2 does not conform with this semiclassical picture. 

A quite appealing feature of these molecular clusters is 
their high point group symmetry, namely the full Octa- 
hedral group O/i = X i (with 48 elements) and the full 
Icosahedral group l/i = I x i (with 120 elements) for the 
cuboctahedron and the icosidodecahedron respectively 
(here i denotes the inversion). This allows for a drastic re- 
duction of the dimensionality of the problem. In order to 
fully exploit all symmetry operations we have employed 
a generalization^ of the standard ED technique^ so as 
to treat higher than onc-dimcnsional Irreducible Repre- 
sentations (IR's) also. With this approach, one is able to 
classify the energy levels according to both Sz and the IR 
of the point group while resolving their full degeneracy. 

The remainder of the article is organized as follows. 
In Sec. |TT] we discuss some of the spectral features of the 
present kagome-like nanomagnets (with particular em- 
phasis on localized magnons) and contrast them to typi- 
cal spectra of unfrustrated AFM's. This is illustrated by 
comparing with a simple 12-site bipartite s = 1/2 clus- 
ter. The investigation of the nature of the M = Ms/3 
plateau is presented in Sec. IIIII This includes both ana- 
lytical and numerical results from high order degenerate 
perturbation theory around the Ising limit, the construc- 
tion of the associated effective QDM's and their extrap- 
olation to the Heisenberg limit. We also discuss the case 
of higher s and the connection to the plateau phase of 
Mo72Fe3o. In Sec. IIVI we demonstrate the presence of 
several low-energy rotational bands in s > 1/2 Heisen- 
berg spectra and reveal their semiclassical origin. The 
analysis is based on a careful comparison to the symme- 
try properties of the semiclassical states and follows the 
basic lines of the seminal works of Bernu et alM'-^ and 
Lecheminant et al^^ on this subject in the context of 
the triangular and kagome AFM. Predictions for the cor- 
responding tower of states are also given for the (s > 1/2) 
icosidodecahedron. Our core idea of the presence of sev- 
eral rotational bands due to the large spatial degeneracy 
of the classical states is also exemplified in Sec. lIVBI bv 
a discussion of the much simpler case of the s > 1/2 
XY model. We leave Sec. |V] for an overview of the ma- 
jor findings of this work. In order for the manuscript 
to be self-contained, we include two appendices. In Ap- 
pendix]^ we summarize the main aspects of the degener- 
ate perturbation expansion around the Ising limit, while 
in Appendix |B] we give the details of the derivation of 
the full symmetry properties of the semiclassical towers 
of states for the Heisenberg and the XY model. 

A special remark is in order here regarding our choice 
of presentation of the spectra. Since we are interested 
in the low-energy excitations in the whole magnetization 




FIG. 2: (Color online) Schematic representation of the un- 
frustrated, bipartite 12-site cluster discussed in Sec. |lTl Its 
symmetry group is Deh = De x i. 



range (these are the most accessible and thus most rele- 
vant as one ramps up an external field at low tempera- 
tures) and in order to best illustrate the central features, 
we have chosen to (except for Fig. [3]) shift the lowest en- 
ergy Eo{Sz) (or Eo{S)) of each Sz (S) sector to zero. This 
guarantees a fine resolution of the low-energy spectra in 
the whole magnetization range. 



II. UNFRUSTRATED VS. KAGOME-LIKE 
AFM'S : GENERAL SPECTRAL FEATURES 

(S=l/2) 

Our main purpose in this section is to present the low- 
energy spectra of the s = 1/2 Heisenberg cuboctahedron 
and icosidodecahedron and to highlight their main fea- 
tures which are common in all frustrated AFM's. For 
comparison, it is expedient to also present the energy 
spectrum of a bipartite unfrustrated magnet. To this 
end, we have chosen the hypothetical 12-site cluster de- 
picted in Fig. [21 The symmetry of this cluster is the 
dihedral group Dg/i = Dg x i which consists of 24 ele- 
ments. Figure [3] shows the low-energy Heisenberg spec- 
trum as a function of S{S+ 1), classified according to the 
12 different IR's of D^h (cf. Ref. 37) shown in the legend. 
The spectrum is typical of unfrustrated AFM's^'^^ of fi- 
nite size N. For instance, we may associate the lowest 
energy band indicated by the dotted line in Fig. [3] with 
the so-called Anderson tower of states^, which is the 
finite-size manifestation of the SU(2) symmetry break- 
ing process occurring in the thermodynamic limit. As 
can be seen in Fig. [31 this tower consists entirely of the 
two one-dimensional representations Alg and B2g of Dg/i, 
which alternate between even and odd S respectively. 
The physics behind this symmetry structure is intimately 
related to the symmetry properties of the semiclassical 
two-sublattice Neel state. For instance, the combina- 
tions Alg ± B2g transform into each other in exactly 
the same way as the two spatial counterparts of the Neel 
state. Above the lowest tower of states of Fig. [31 there 
exists a finite excitation gap followed by a quasi contin- 
uum of higher excitations. All these features are typical 
of unfrustrated AFM's. 

In contrast, frustrated AFM's show very different low- 
energy features as exemplified by the s = 1/2 spectra 
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FIG. 3: (Color online) Low-energy spectrum of the 12-site 
bipartite AFM shown in Fig. [2] as a function of S{S + 1) and 
classified according to IR's of the Dgh group. The dotted 
line denotes the Anderson tower of states which embodies the 
finite-size features of the spatial and SU(2) broken Neel state 
in the thermodynamic limit. 



shown in the two lower panels of Fig.Ufor the two molec- 
ular magnets of the present study. For comparison, the 
upper panel shows the low-energy portion of Fig. [3] in 
terms of 5*2. The contrast between the two types of spec- 
tra is more than evident (note that both of the upper 
two panels correspond to 12-site clusters and are shown 
in the same energy scale). The most striking feature 
emerging in frustrated AFM's is the absence of a clear 
energy scale separating a lowest band from higher lying 
excitations. Instead, a large "bulk" of low-energy exci- 
tations is manifested in the whole range of Sz forming 
a quasi-continuum. This is a central feature that holds 
also for higher s (cf. Sec. IIV[) and stems from the highly 
frustrated exchange interactions in these clusters. The 
nature of these excitations for s — 1/2 is not completely 
understood^ but, as we are going to show in Sec. IIV[ a 
qualitative understanding can be obtained for s > 1/2 
based on the large classical degeneracy of spin configura- 
tions which remains dominant in the semiclassical regime. 
In particular, the broad INS response reported in Ref.|28l 
for Mo72Fe3o is naturally accounted for by the results of 
this analysis. 

Let us now describe shortly some special spectral fea- 
tures and their origin. The ground state energies Eq{Sz) 
for the two nanomagnets for s = 1/2 are given in Table [H 
These energies determine the zero-temperature magneti- 
zation processes shown in Fig. [51 For the excitations 
above the ground state, three regimes of special interest 
can be highlighted (see shaded areas in Fig. 2]): (i) the 
singlet excitations below the lowest triplet^- which are 
given in Table [III and amount to 7 for the cuboctahedron 
and 80 for the icosidodecahedron^, (ii) the existence of 
degenerate localized magnons below saturation and (iii) 
the presence of a number of well isolated low-lying states 
right above the M = Ms/3 {Sz = 5) plateau phase of the 
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FIG. 4: (Color online) Low-energy spectra (shifted as de- 
scribed in the text) of the s — 1/2 Heisenberg model on the 
12-site unfrustrated magnet shown in Fig. [2l (top), on the 
cuboctahedron (middle) and the icosidodecahedron (bottom). 
Three special features are highlighted by the corresponding 
shaded areas in the two lower panels: (i) the low-lying sin- 
glets below the first triplet in the Sz = sectors, (ii) the 
existence of localized magnons highlighted in the sectors be- 
low saturation, and (iii) the lowest 36 Ising-like configurations 
(cf.Sec. lIlIB2|) above the plateau Sz = 5 sector of the icosido- 
decahedron case (lowest panel). The large energy gap between 
these configurations and higher excitations is indicated by the 
arrow. 



TABLE I: Lowest energies Eo{Sz) of each Sz sector and the 
corresponding degeneracies for the s = 1/2 cuboctahedron (a) 
and icosidodecahedron (b). 
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TABLE II: (a) Energies of the seven lowest singlets of the 
s = 1/2 Heisenberg cuboctahedron lying below the first triplet 
E = -5.06220685J (T2u) state, together with their Oh sym- 
metry classification and their degeneracy, (b) Energies of the 
80 lowest singlets of the s — 1/2 Heisenberg icosidodecahe- 
dron lying below the first triplet E = -13.01640033 J (Ag), 
together with their 1^ classification. 
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(b) s = 1/2 icosidodecahedron 
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s — 1/2 icosidodecahedron. The latter will be analyzed 
in detail in Sec. IIIII 
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FIG. 5; (Color online) Zero-temperature magnetization 
curves for the Heisenberg cuboctahedron (s — 1/2, 3/2, 5/2) 
and the s — 1/2 icosidodecahedron (see also Refs. Ll5,18). 
The saturation magnetization and field values are given by 
Ms — Nss{gjj,B) (where Ns is the number of sites) and 
Hs = 6sJ/(g^s). 



The concept of localized magnons has been largely 
discussed in the context of highly frustrated bulk 
ArM's^i^^i^i^iiiii^i^. For the present clusters it has 
been discussed by Schnack et ali^^^. We shortly re- 
visit this issue here in the light of our symmetry re- 
solved method. Quite generally, the eigenvalues of H 
in the one magnon {M — Ms — 1) subspace are equal 
(apart from an overall constant energy shift) to the eigen- 
values of the adjacency matrix {J^^u) times the spin 
s^. Our decomposition of the respective subspaces for 
the cuboctahedron and the icosidodecahedron in terms 
of IR's of the full Oh and 1^ groups are (in order 
of increasing energy): (Eg®T2u)eT2g©TlueAlg, and 
(Hg©Hu)©T2u©Fg©Fu©Hg©Tlu©Ag respectively, and 
they are compatible to the ones given in Ref . |4J in terms 
of IR's of the and I subgroups. 

For the cuboctahedron, the lowest one-magnon level is 
5-fold degenerate (Eg©T2u), see Fig. gl^ middle) . These 
correspond to localized, non-interacting magnon states. 
The smallest loops that can host such magnons are the 




FIG. 6: (Color online) The minimal loops hosting 
the localized k — ir magnons on the topology of 
the cuboctahedron (a) and the icosidodecahedron (b). 
They read \ip) = 5 (s^ — Sj" + s^ — sj) |0), and \ip) = 
^75 {si — sj + . . . — Sg) |0) respectively, where |0) is the fer- 
romagnetic vacuum. These are exact eigenstates of Eq. ^ 
owing to a cancellation of interaction terms resulting from the 
special corner-sharing triangles topology"^- . In both cases, the 
energy hujm required to excite these localized states measured 
from the ferromagnetic (FM) vacuum |0) (-Bo = NbJs'^, where 
Nb is the number of bonds) equals huim = — 6sJ, and is inde- 
pendent of the length of the loops. 



The existence of localized, non-interacting niagnon 
states results in a magnetization jump of AS 2: > 1, since 
the lowest energies at the corresponding Sz sectors scale 
linearly with the number of magnons, and thus cross 
each other at the same (saturation) field. We remark 
here that all features related to the existence of local- 
ized magnons (symmetry decomposition, degeneracy, and 
the magnetization jump in absolute units) do not depend 
on the value of s (see e.g. Fig. [14] below). Finally, the 
fact that the number of independent magnons is larger 
in the icosidodecahedron than the cuboctahedron case 
is clearly related to their size. In extended frustrated 
AFM's, this number grows exponentially with system size 
but depends in a non-trivial way on the topology of the 
system and is connected to the question of linear inde- 
pendence^!^. The extensive degeneracy gives rise to a 
macroscopic magnetization jump at the saturation field 
and a large magnetocaloric effect (see e.g. Ref. l4lr ). A 
study of the latter on the present clusters can be found 
in Ref. M. 



square faces depicted in Fig. [6ja). The 5- fold degener- 
acy is due to the fact that there are 6 different square 
faces on this cluster, but not all magnons are indepen- 
dent: The sum of all 6 square magnons taken with oppo- 
site phases in neighboring squares vanishes. The lowest 
energy level of the Sz = 4 two-magnon space is 3-fold 
degenerate (Ag®Eg), and corresponds to the 3 different 
ways of placing two non-interacting magnon excitations 
(there are three different pairs of opposite squares). Plac- 
ing one more magnon gives an interaction energy cost and 
a non-degenerate Sz — 3 lowest level. We should remark 
here that niagnon states "living" on the hexagonal equa- 
tors of the cluster are also exact eigenstates, but each of 
these can be easily expressed as a linear combination of 
surrounding square magnons. The lowest level degenera- 
cies for the one-magnon and the two-magnon space are in 
agreement with the values of 7V/3-t-l and N'^/18-N/2+l 
respectively with A^ = 12 which are derived in Ref. |4Q for 
the kagome lattice (for which the hexagonal loops are the 
most local and natural ones for the description of the lo- 
calized magnons). 

For the icosidodecahedron, the lowest one-magnon 
level is 10- fold degenerate (Hu®Hg). Here, the small- 
est loops that can host such localized states are the 
octagons surrounding a given vertex and depicted in 
Fig. IH^b). The 10- fold degeneracy can be attributed to 
a non-trivial linear dependence among the 30 different 
octagonal magnons on this cluster. The lowest energy 
level of the two-magnon manifold is 25-fold degenerate 
and decomposes into Ag © Au © Fg © Fu ® 2Hg © Hu. 
Hence, there exist 25 ways of placing two mutually non- 
interacting magnons. Similarly, the lowest energy of 
the three-magnon space is two-fold degenerate (Ag©Au), 
whereas that of the S'^ = 11 sector is non-degenerate, 
signifying that it is not possible to have four magnons 
without an interaction energy cost. 



III. M = Ms/3 PLATEAU PHASE 

In this section, we focus on the nature of the exci- 
tations above the M = Ms/3 plateau. There are two 
major reasons for paying special attention to this partic- 
ular plateau among the remaining ones which are present 
anyway in our finite-size clusters (cf.Fig. [5]). The first is 
of practical interest and is related to the experimental 
manifest at ionil of this particular phase in the s — 5/2 
Mo72Fe3o cluster. Besides, as shown in the upper panel of 
Fig. [51 the M = Ms/3 plateau seems to be the most sta- 
ble and survives at finite s > 1/2 (the staircase s = 1/2 
magnetization process eventually turns into the expected 
(classical) linear behavior-il for very large s) . The second 
reason is that the M = Ms/3 plateau phase is a generic 
feature of frustration and is known to survive in the ther- 
modynamic limit for some bulk AFM's (cf. Ref. [21). 

In Sec. IIII Al we present and analyze a striking fea- 
ture of the excitations above the plateau phase of the 
s = 1/2 icosidodecahedron and show how it can be ob- 
served experimentally in thermodynamic measurements. 
In Sec. IIIIBI we present our derivation of an effective 
Quantum Dimer Model for the plateau and reveal the 
major role of the topology and the spin s. 



A. Thermodynamics 

Our Exact Diagonalizations for the s ~ 1/2 icosido- 
decahedron shown in the lowest panel of Fig. [3| reveal a 
striking feature at the Af = Ms/3 sector: The low-lying 
excitation spectrum immediately above the plateau con- 
sists of a group of 36 states (their energies are given in 
Table Hill) which are well isolated from higher excitations 
by a gap of order 0.2 J, an order of magnitude larger than 
the excitations {^ O.OIJ) within this manifold. We ar- 



TABLE III; Energies (in units of J) of the 36 lowest S^ = 5 
states of the s — 1/2 Heisenberg icosidodecahedron together 
with their \h classification. The lowest excitation above this 
manifold lies at _B = -10.04843786 (Hg). 



Energy [J] 


IR(deg) 


Energy [J] 


IR(deg) 


-10.30278978 


Hu (5) 


-10.26904953 


Hu (5) 


-10.29875816 


Ag(l) 


-10.26657194 


Hg(5) 


-10.29837409 


Hg(5) 


-10.25765943 


Hg(5) 


-10.29057364 


Fg(4) 


-10.25604215 


Ag(l) 


-10.28622445 


Fu (4) 
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FIG. 7: (Color online) Temperature dependence of the mag- 
netic entropy Sm and specific heat Cm (in units of Boltz- 
mann's constant ks) content of the lowest 36 plateau states 
of s = 1/2 Heisenberg icosidodecahedron. Here the magnetic 
field corresponds to the center of the M = Ms/3 plateau. The 
entropy starts off from the value ln5 ~ 1.6094 corresponding 
to the lowest 5-fold degenerate Hu state. The dashed line 
denotes the entropy value of In 36 ~ 3.5835 corresponding to 
the full 36-fold low-energy subspace. 



gue below that this peculiar feature must be related to 
the special topology of the icosidodecahedron (it does not 
appear for the cuboctahedron). Given this behavior, it 
is expedient to consider the low-temperature dependence 
of the magnetic specific heat and entropy content of the 
lowest 36 plateau states. These quantities are shown in 
Fig. [7] for the s = 1/2 case and for giisH ~ 1.1326J 
(where g is the electronic spectroscopic factor and /is 
the Bohr magneton) which corresponds to the center of 
the M = Ms/3 plateau. At this field value, the lowest 
excitations of the adjacent Sz — 4: and Sz = Q sectors lie 
approximately 0.2J above the ground manifold (higher 
Sz = 5 states lie ~ 0.25 J above). Hence, the temperature 
behavior shown in Fig. [7] must be valid at fc^T < O.IJ. 
In this temperature regime, the entropy content of the 
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FIG. 8: (Color online) A local view of one of the "uud" Ising 
configurations (with two spins pointing up and one down in 
each triangle) and the mapping between vertices of the cuboc- 
tahedron (solid line in (a)) and edges of its dual cluster, the 
cube (dashed lines). In (a), all spins point up except the ones 
at vertices 1 and 3 (designated by the black dots) which point 
down. By mapping each down spin in (a) to a dimer on the 
corresponding edge of the cube we obtain the dimer plaquette 
in (b). As discussed in Sec. HUB II such square loops with al- 
ternating up-down spins have the minimum even length and 
thus govern the lowest order kinetic processes driven by Hxy ■ 
In (a), these read t sfs2 s'^sj with t ex a*" and map to the 
dimer plaquette fiip of Eq. pi|l . The lowest order diagonal 
processes are also confined on these square loops and scale as 



lowest Sz = 5 states is already saturated to its full value 
of In 36 ~ 3.5835 (dashed line) which amount to a siz- 
able fraction of about 17% of the fuU 30 x In 2 ~ 20.7944 
magnetic entropy of the cluster. The fine details of Fig. [7| 
can be associated to the actual splitting between the 36 
states. For instance, the double-peak form of the specific 
heat stems from the small separation of the first 19 from 
the remaining 17 states (cf.Table lllll and lowest panel of 
Fig. 2]). Note also that the entropy starts off from the 
value ln5 ~ 1.6094 corresponding to the lowest 5-fold 
degenerate Hu state (cf.Table lHip . 



Unfortunately the plateau regime of M072V30 cannot 
be reached experimentally due to the large exchange 
value of J/ks — 250 K. We shall argue below, based on 
the results of our effective QDM, that a similar structure 
must exist above the M = Ms/ 3 plateau of the s = 5/2 
Mo72Fe3o cluster. In particular, we shall argue that (i) 
the lowest 36 states are split into two almost degener- 
ate levels of 30 and 6 states respectively with the former 
being lowest in energy, and (ii) that a sizable gap be- 
tween these 36 states and higher excitations must proba- 
bly survive as well. The corresponding specific heat peak 
can be verified by thermodynamic measurements on the 
Mo72Fe3o cluster at the plateau regime of H ~ 5.9 Tesla. 
Despite the very small exchange value (J ~ 1.57 K) of 
Mo72Fe3o (one presumably needs to reach ultra low tem- 
peratures, T < 200 mK) one may still confirm our picture 
by an assessment of the missing entropy-^5:-^. 




where 
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FIG. 9: (Color online) A local view of one of the "uud" Ising 
configurations (with two spins pointing up and one down in 
each triangle) and the mapping between vertices of the icosi- 
dodecahedron (solid line in (a)) and edges of its dual cluster, 
the dodecahedron (dashed lines). In (a), all spins point up ex- 
cept the ones at vertices 2, 4, 6 and 8 which point down. By 
mapping each down spin in (a) to a dimer on the correspond- 
ing edge of the dodecahedron one obtains the dimer plaquette 
in (b). As discussed in Sec. HUB 21 such octagonal loops with 
alternating up-down spins have the minimum even length and 
thus govern the lowest order kinetic processes driven by Tixy 
In (a), these read t s^ s^s^ ■ ■ ■ ^t with t oc o?" and map to 
the dimer plaquette flips of Eq. (|16|l . On the other hand, 
the lowest order diagonal processes are confined in a single 
pentagon and scale as u ex a^. 



B. Effective QDM 

In what follows we present a complementary picture 
for the nature of the excitations above the plateau phase. 
This picture reveals the central role of the topology and 
the intrinsic spin s and will emerge from the derivation 
of an effective quantum dimer model (QDM) in the spirit 
of Refs. [gl llOlilll . The main idea is to start from the de- 
generate M/Ms — 1/3 ground state configurations of the 
Ising limit and establish an adiabatic connection to the 
low-lying excitations of the Heisenberg point by employ- 
ing a perturbative expansion in the anisotropy parameter 
CK = Jxy/Jz- The resulting effective Hamiltonian can be 
cast into the form of a QDM on the dual clusters as ex- 
emplified in Figs. [5] and [5] for the cuboctahedron and the 
icosidodecahedron respectively. For a general spin s, the 
M = Ms/3 ground state manifold of the Ising Hamilto- 
nian Tiz spans all configurations with two spins having 
m — s and one with m — ~s in each triangle. Each 
one of these "uud" states on the cuboctahedron and the 
icosidodecahedron is in one-to-one correspondence to a 
closed-packed dimer covering on their dual clusters, the 
cube and the dodecahedron respectively. 



1. Cuboctahedron 

We start with the "uud" GS's of the Ising cuboctahe- 
dron. There are nine such states since this is the number 
of different dimer coverings on the cube. This manifold, 
henceforth Pinvir]^, decomposes into two invariant (under 

Puud ^it^ 6 ^^"i "^ states 



Oh) "uud" famines pI.'^^^ and 
respectively: 



>(6) 



(3) 



Pu.^ = Alg ©Eg © T2u, P];^^ = Alg © Eg . (6) 

Each of the two families contains states with a fixed num- 
ber Uc (2 and 4 respectively) of square plaquettes of the 
type of Fig. IllT c) (this number remains invariant under 
the operations Oh of the cluster). 

The 9 "uud" states are highlighted in Fig. [TO] which 
shows our symmetry-resolved ED results for s — 1/2, 
3/2, and 5/2 (at their M — Ms/3 sector) as a function 
of the anisotropy parameter a. The energies are given 
in units of JzCt^ and JzCt"^ for s = 1/2 and s = 3/2, 5/2 
respectively, which are the leading orders of the energy 
splitting due to Hxy (see below) . Figure [TO] shows that 
the Heisenberg states which are adiabatically connected 
to the lowest Ising manifold are not the lowest excitations 
for s — 1/2 while this is clearly the case for higher spins 
and, as we show below, for the s = 1/2 icosidodecahedron 
as well. 

We shall try now to understand some of the features 
of Fig. [TO] in more detail by considering the lowest order 
effect of T-Cxy in splitting the Ising nine-fold degenerate 
manifold, as a function of spin and a. We follow the gen- 
eral guidelines and considerations of the Appendix]^] We 
distinguish between diagonal and off-diagonal processes 
depending on whether the initial state is finally recovered 
or not. The former come from the smallest closed paths 
on the molecule (beyond triangles) , which in the present 
case are the square loops (see left panel of Fig. [1]) and the 
corresponding amplitudes scale with the fourth power of 
a for all s. Now, there are only three possible configu- 
rations on a square which respect the "uud" constraint 
and these are depicted in Fig. [TTT a). (b) and (c). Each 
one carries a certain diagonal energy, say Ca, £& and tc- 
We have calculated these energies as a function of s using 
Eq. (|A2p and by enumerating all relevant processes. The 
results are 



ea(s) 
efc(s) 

ec(s) 



= 



2(4s-l)2 



Jza , 



/2 



(4s- l)2(2s-l) 



J.a' 



(7) 



uud 



,(6) 



,(3) 



(5) 



The Kronecker symbol Ss^i/2 appears in eds) because 
some of the diagonal processes relevant for s > 1/2 are 
not present for s — 1/2 since they involve intermediate 
states which do not belong to the Ising manifold. ^^ Now, 
to any given Ising configuration i there corresponds an 
associated potential energy equal to E^ — n\ea + n^e^ -f 
n\ec^ where n^, n^, n* are the number of squares in the 
states a, b and c, respectively in i. On the other hand, 
we must satisfy two global conditions, one for the total 
number of squares Ns = 6 = nj, -|- n^ -|- rt* , and another 
for the total number of down spins N^ = 4 = n\/2 + n],. 
This leaves us with one independent, non-global variable, 
say n^, in terms of which one can express E^ . Omitting 
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FIG. 10: (Color online) Lowest energy spectrum of the XXZ 
model on the cuboctahedron as a function of a = Jxy/Jz, 
for s = 1/2 (top), 3/2 (middle), and s = 5/2 (bottom). The 
energies are given in units of JzCt^ (top) and J^a^ (middle and 
bottom), which correspond to the leading orders of the energy 
splitting due to Tixy (see text). For s = 1/2, the dominant 
off-diagonal processes split completely the 9-fold degeneracy 
of the Ising point. For higher spins, the diagonal processes 
dominate and give rise to a splitting of 2v{s,a) between the 
two "uud" families of Eq. ^ as predicted from perturbation 
theory (see text) . Filled squares denote the eigenvalues of the 
effective QDM of the corresponding leading term for each s. 
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FIG. 11: The three possible local (on the square loops) views 
of the "uud" configurations of the cuboctahedron. Filled cir- 
cles denote spins with m = — s, all remaining vertices have 
m = s. 



a global energy term 8eb — 2ea, we find E^ = v{s,a)nl 
with 



u(s, a) 



2eh + Cc 



s^{2sds,i/2 - 1) 
(4s-l)2(2s-i; 



J,a^ , (8) 



which is positive for s — 1/2 and negative otherwise. The 
corresponding (lowest order) effective diagonal Haniilto- 
nian reads 



v'.^^=v{s,a)Y: 



(9) 



where the sum runs over all six square plaquettes (with 
both orientations) of the cluster. We may easily check 
that among all nine possible dimer coverings of the cube, 
three of them have tt-c = 4 and thus Vgg = 4v{s,a), while 

the remaining six have Uc = 2 and thus V^J = 2v{s, a). 
These correspond to the two families of "uud" states 
mentioned above, see Eq. ([5]). Hence, the eigenvalues of 



y 



^g form a pair of a 6-fold and a 3-fold degenerate lev- 
els with an energy splitting of 2v{s, a) between them. In 
particular, this splitting amounts to -^JzCt^ for s = 3/2 
and j^Jza"^ for s = 5/2. 

We now consider off-diagonal processes. To lowest or- 
der in a, these are confined to the maximally flippable 
even-length loops of the cluster. These loops are the al- 
ternating spin up-down configurations already shown in 
Fig. [H The corresponding flipping amplitude t scales as 
a^*. Their explicit values for several s are provided in 
Table En] (with L = 4) of Appendix E] Thus the leading 
kinetic effect is described in the dimer representation by 
the term 



T-(4s) 
-'eff 



^(^'")E(|D>(D 



(10) 



where i(s, a) can be calculated explicitly using 
Eq. (|A2[) and enumerating all different processes. 
Several representative values are provided in Ta- 
ble IVIII The eigenvalues of T^g* in units of t 
are: -2\/2, -\/2, -\/2, 0,0, 0, \/2, V^, 2\/2. These cor- 
respond to a complete splitting of the different IR's of 
each of the two families of Eq. ^. 

According to the above, our effective quantum dimer 
model should generally include both kinetic and potential 
terms. To leading order, this model reads 



Heff = V, 



(4) 
eff 



T-(4s) 
-'off 



(11) 



10 



For s = 1/2 we have t = —J^a^ 



Jza' 



and the 



dynamics is mainly governed by kinetic processes which 
spht completely (cf. Fig. [TOTa)) the nine "uud" states. 



For s = 1 we have t — —J, a 



-Jztt /9 and thus 



both diagonal and ofF-diagonal processes are equally im- 
portant. For s > 1 the diagonal processes dominate and 
give rise to a splitting of 2u(s, a) between P^^tj and P^^j 

of Eq. ([5]). In particular, since u < 0, the states of P^^j 
will be favored because they have a larger number (four) 
of the plaquettes of Fig.[8l[c). All these features are nicely 
demonstrated in Fig. [TU] where we compare our ED re- 
sults for the XXZ model at small a with the leading 
order eigenvalues of Eq. pT|) which are shown as (red) 
filled squares. 



2. Icosidodecahedron 

We turn now to the corresponding plateau phase of the 
Heisenberg icosidodecahedron and follow a similar anal- 
ysis as above. Here the lowest Ising manifold, henceforth 
■^uud' consists of 36 "uud" states which are in one-to- 
one correspondence with the 36 dimer coverings of the 
dodecahedron. We find that this manifold decomposes 
into two invariant (under l/j) families R-uu^j and R[j^(j of 
30 and 6 states respectively as 



R 



uud 



R 



(30) 
uud 



R 



(6) 
uud 



where 



r|^3^o) ^ Ag©Au®Fg©Fu®2Hg®2Hu 



R 



(6) 
uud 



= AgeHg 



(12) 



(13) 



Each of these families contains states with a fixed number 
TT-c (8 and 10 respectively) of pentagonal plaquettes of 
the type of Fig. I12r c) (this number remains invariant 
under the symmetry operations U of the cluster). We 
should note in particular that each of the 6 states of 



R 



(6) 
uud 



contain two pentagons with all spins pointing up 



(i.e. have Ua — 2, cf.Fig.[T2ja)). 

A simple inspection of Fig. [T5f a). which shows the 
s — 1/2 lowest energy spectrum of the XXZ model, re- 
veals that the lowest 36 Heisenberg states trace back to 
the ground state "uud" manifold of the Ising point. A 
striking difference to the s = 1/2 cuboctahedron case 
studied above is that as these lowest Ising states "evolve" 
toward their low- lying Heisenberg counterparts, they re- 
main always well separated from the higher energy states 
in the Sz = 5 subspace. We now give a complementary 
picture, which is valid at least for small a, by considering 
the lowest order processes driven by Tixy and by deriving 
the corresponding effective QDM on the dodecahedron. 
We begin with the lowest order off-diagonal processes. 
As above, these stem from maximally flippable loop con- 
figurations of the smallest possible even length L. Such 
a loop configuration that respects the local "uud" con- 
straint is the octagonal loop with alternating up-down 




FIG. 12: The three possible local (on the pentagonal loops) 
views of the "uud" configurations on the icosidodecahedron. 
The filled circles denote spins with m = ~s, all remaining 
vertices have m = s. 



spins depicted in Fig. [Sfa) which in turn maps to the flip- 
pable plaquette of the dodecahedron shown in Fig. [D^b) . 
For spin s = 1/2 then, the lowest off-diagonal term in 
HeB is of fourth order in Hxy Each "up-down" loop of 
the type shown in Fig. [Ifa) is amenable to a kinetic pro- 



cess of the form t s^ st s 



From our calculations, 

„4 



1 '='2 ''3 • ■ ■ ■'8 

shown in Table rvTTl we obtain t — ~2.bJzCt^- In fifth 
order, we find two types of kinetic processes. The first is 
similar to the above but now involves loops of length 10 
such as the equators of the molecule. The second type 
is less obvious, and invokes again the octagonal loops of 
Fig.lHJa) and any one of the neighboring spin sites. Since 
these loops map to exactly the same flippable dimer pla- 
quette of Fig. [9Ub) their effect is to merely renormalize 
the fourth-order amplitude t. In fact, these terms re- 
sult in an overall decrease of |t| since they carry an extra 
negative sign. 

On the other hand, the lowest order diagonal processes 
must be confined to the smallest closed path which in 
the present case are the pentagons of the molecule (see 
right panel of Fig. [T]). The three possible types of con- 
figurations that respect the "uud" constraint around a 
pentagon are depicted in Fig. [T^] and are designated by 
a, h and c. Each one carries a certain diagonal energy, 
say Ca, Cfc and Cc. These are calculated by enumerating 
all relevant processes and using Eq. (|A2|) of Appendix [XI 
They are explicitly given by 



= 



Cfc 



Cc = 



8(4s-l)2 
6s4 



Jza^ , 



(4s-l)(8s-3): 



rJ.a^ 



(14) 



To any given Ising configuration i there corresponds an 
associated energy equal to E^ = n^^ea + nleb + nlec, where 
n^, nl, n^ are the number of pentagons in the states a, b 
and c respectively in i. On the other hand, we must again 
satisfy two global conditions, one for the total number of 

and one for the total 
. This leaves 

us with one independent, non-global variable, which we 
choose to be n* . Omitting a global energy 20e6 — 8ea, we 
find E^ = v{s,a)nl., with 



pentagons Np = 12 = nl^+nl-\ 
number of down spins Nj = 10 = nl/2 - 



v{s, a) 



-2eb 



s3(32s2 + 24s-9) 

4(4s-l)2(8s-3)2 



J.a', (15) 
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FIG. 13: (Color online) (a) Lowest eigenvalues of the s — 1/2 
XXZ model on the icosidodecahedron in the Sz ~ 5 sector. 
Interpolation between the a = Ising and the a — 1 Heisen- 
berg point. The large energy separation between the lowest 
36 states and higher excitations is clearly evident, (b) Con- 
vergence of the lowest 36 eigenvalues (in units of J^a'^ which 
is the leading order) toward the eigenvalues (filled squares) of 



(4) 



the effective dimer Hamiltonian 7^^ 
in the text. 



of Eq. (|16|l as described 



which is positive for all s. This means that v(s, a) favors 
configurations with the minimum number of the pentag- 
onal states of Fig. flW c). Since all 30 states of R^^j have 

ric = 8 while the 6 states of R^u^ have ric — 10, the 
former family will be lower in energy by a splitting of 
2v{s, a). Furthermore, it is clear that diagonal processes 
do not give rise to a splitting within the two families. 

s = 1/2 case. — Given all the above, the effective QDM 
for the plateau phase of the s = 1/2 Heisenberg icosi- 
dodecahedron is, at lowest order, dominated by kinetic, 
off-diagonal processes of the form 



n, 



off ■ 



r 



(4) 



off 



= t 



z- \ 



\\// 



+ h.c. 



(16) 



where the sum runs over all octagonal plaquettes and 
t — —2.bJzOi'^ ■ The corresponding 36 x 36 Hamiltonian 
matrix can be constructed and solved numerically for its 
eigenvalues. These are provided in Table IIVI in units of 



TABLE IV: Eigenvalues (in units of t) of Ticff given in 
Eq. (|16p . together with their multiplicities. 



Energy [f] 


deg. 


Energy [f] 


deg. 


-4 


1 


-0.694593 


5 


-3.06418 


5 


1 


4 


-2.89898 


1 


2 


5 


-2 


5 


3.75877 


5 


-1 


4 


6.89898 


1 



t. They are also shown as filled squares in Fig.fTWb). In 
the same figure, we show the lowest 36 eigenvalues (di- 
vided by Jz) of Ti' in units of a^. The clear convergence 
for small a toward the eigenvalues of Ticft confirms the 
validity of our perturbative calculations. Moreover, the 
fact that the convergence is linear confirms that the next 
processes contributing to TLcs come in fifth order. In 
particular, the clear decrease of the bandwidth with a is 
in agreement with our previous assertion that the fourth 
order amplitude t of Eq. p6p gets renormalized from the 
fifth order octagonal kinetic processes mentioned above. 
The latter seem to dominate over the corresponding fifth 
order off-diagonal decagonal loop processes and the fifth 
order diagonal ones. Looking at Fig.fT^Ta) one notes that 
this may be even more general: To all orders in TLxy^ 
there seems to be a mere renormalization of the band- 
width without drastically altering the relative amplitudes 
of the eigenvalues of T^^ ' . This suggests a dominance of 
the most local (octagonal) kinetic processes renormalized 
from all orders in Tixy 

s > 1/2 and relevance to Moj2Fe^Q. — As mentioned 
above, diagonal processes first appear in fifth order irre- 
spective of s. On the other hand, the lowest off-diagonal 
process on a loop of L (even) sites appears in order Ls. 
For instance, the octagonal loops discussed above give 
processes at order 8s (and higher), whereas decagonal 
(e.g. the equatorial) loops contribute in order IDs (and 
higher). Hence for s — 1 diagonal and off-diagonal pro- 
cesses are equally important, while for s > 1 the physics 
will be completely dominated by diagonal processes of 
fifth order in a: 



n, 



off ■ 



v. 



(5) _ 



v{s, a) 



i^W 



(17) 



where the sum runs over all pentagonal plaquettes (with 
all possible orientations) and v{s, a) is given by Eq. (|15p . 
As explained above these processes result in a diago- 
nal sphtting of 2u(s, a) between R„„j and R^Jd with 
the former family being lowest in energy. For s = 5/2 
in particular, which is relevant for the plateau phase of 
Mo72Fe3r/'', this splitting amount to the quite small value 
of 2v ~ 0.0838 JzOi^ ■ Since, at the same time, the exci- 
tations out of the Ising manifold are expected to remain 
gapped for the finite value of s = 5/2, we speculate that 
the plateau phase of Mo72Fe3o must show a characteristic 
low-temperature thermodynamic signal which is qualita- 
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tively similar to that of Fig. [7] (up to fine details related 
to the 30 to 6 diagonal splitting). Here, in particular, 
the entropy content of the renormalized "uud" manifold 
amounts to approximately 7% of the full magnetic en- 
tropy. This calls for low-temperature specific heat mea- 
surements on Mo72Fe3o at the plateau phase as explained 
in Sec. IIII Al More generally, it is exciting that the no- 
tion of a Quantum Dimer Model finds a realization in 
the low-energy M = Ms/3 plateau physics of finite-size 
magnetic clusters like Mo72Fe3o or M072V30. 



IV. HEISENBERG SPECTRA FOR s > 1/2 

The major focus of this section is on Heisenberg spec- 
tra with s > 1/2. Our interest in this regard is mainly 
motivated by the INS experiments reported by Garlea et 
al^ on Mo72Fe3o. The main finding of these experiments 
is a very broad response which manifests in a wide range 
of fields. Previous theories which are based either on the 
excitations of the rotational band mode l^^'^^ or on spin 
wave calculation o^°'^^ , could not account for the observed 
behavior since they predict only a small number of dis- 
crete excitation lines at low temperatures. Although the 
diagonalization of the s > 1/2 icosidodecahedron is not 
feasible (at low magnetizations) with current computa- 
tional power, an immediate interpretation of this behav- 
ior can be deduced from a study of the cuboctahedron. 
Exact diagonalization spectra of this cluster for s = 3/2 
and s = 5/2 are shown in the two lower panels of Fig. [Til 
For comparison we have also included the s = 1/2 spec- 
trum (this was shown before in the middle panel of Fig. [4] 
in terms of Sz instead of the total spin S). Here, in con- 
trast to unfrustrated clusters (cf. upper panel of Fig.H]), 
there does not exist any well isolated and thus clearly 
identified low-energy tower of states or rotational band. 
Instead, a "bulk" of very dense excitations are present in 
the full magnetization range. This is a generic feature of 
highly frustrated systems which manifests irrespective of 
s and thus must be also present in the s — 5/2 Mo72Fe3o 
cluster (similarly to its s = 1/2 analogue of Fig. [5]). 

The origin of these dense excitations for s > 1/2 can 
be readily suggested by the following striking observa- 
tion in Fig. [131 The spectra consist, up to a relatively 
large energy cut-off, entirely of the representations Alg, 
A2g, Eg, and T2u. The main message in the following 
is that this peculiar spatial symmetry pattern as well as 
the combined spatial -|-spin pattern (i.e., the appearance 
of specific sets of spatial IR's in each S sector) are char- 
acteristic fingerprints of the 3-sublattice Heisenberg clas- 
sical GS's. For instance, as we explain below, these four 
representations are exactly the ones that appear in the 
symmetry decomposition of the coplanar classical GS's 
(cf. Eq. p^ below) . The dense excitation features of the 
lower two panels of Fig. [M] can be thereby accounted for 
by the large spatial degeneracy of these configurations, a 
fact whose importance does not seem to have been rec- 
ognized in the past. In principle, each of these states 



gives rise to a distinct "tower of states" or "rotational 
band" and they all appear together at low energies al- 
beit split by quantum fluctuations. So the large discrete 
degeneracy has a direct impact on the low-energy spec- 
trum, and we believe this large number of levels is at the 
very heart of the broad INS response reported in Ref. |28| . 
By contrast, the absence of a clear symmetry pattern in 
the s = 1/2 spectra (cf. upper panel of Fig. [H|) shows 
that the associated low-lying excitations are of different 
origin. 

Before analyzing further our numerical results it is use- 
ful to recall (cf. Subsec. IIV Al below) what is known 
about the classical GS's of the infinite kagome lattice 
in zero and finite field and discuss what carries over in 
the present clusters. In particular, we give the explicit 
spatial degeneracy of the coplanar GS's and a short sum- 
mary of their symmetry properties. The latter have been 
derived independently by employing a group theoretical 
analysis, the details of which have been relegated to Ap- 
pendix [Bj Our semiclassical interpretation for the origin 
of the dense excitations of the above s > 1/2 spectra 
will be further corroborated by a closer comparison of 
the symmetry pattern of the spectra with the combined 
spatial -l-spin symmetry of the coplanar GS's also derived 
in Appendix |B] In Subsec. IIVBI we discuss the simpler 
case of the s > 1/2 XY model which exemplifies very 
evidently the main idea of this section, i.e., the simulta- 
neous presence of several lowest towers of states due to 
the discrete classical degeneracy. 



A. Classical GS's and large spatial degeneracy 

Let us first consider the ground state configurations 
of the classical Heisenberg model in the infinite kagome 
system and the present clusters. The corner-sharing tri- 
angles structure makes the discussion rather simple. The 
classical Hamiltonian can be rewritten in the suggestive 
form (in units of gfiB ~ 1) 



K 



HB 

classical 



J 



= - ^(Sa-H/2J)^ 



(18) 



where Sa denotes the total spin on a triangle A. In this 
form it is straightforward to see that all configurations 
with Sa = H/2J on each triangle are GS's. It is useful 
to examine the zero- field case first. 

Zero-field case. — Here, the classical constraint Sa = 
amounts to a simple 120° configuration of the three spins, 
which is depicted in Fig. llSf a) . An important point here 
and in the following is to determine how many such GS's 
exist. For the kagome lattice it is well known that the 
ground state manifold consists of both coplanar and non- 
coplanar configurations in zero field. The coplanar GS's 
are extensively degenerate as can be shown by a map- 
ping onto vertex three-colorings of the kagome lattice or 
equivalently onto bond three-colorings of the Honeycomb 
lattice^i^. On the other hand the non-coplanar GS's 
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FIG. 14: (Color online) Low-energy spectra of the AFM 
Heisenberg cuboctahedron for s = 1/2 (top), 3/2 (middle) 
and 5/2 (bottom) in terms of the total spin S. All energies 
are measured from their corresponding Eq{S). 



can be generated from the coplanar ones by the follow- 
ing recipe: Identify a loop of alternating spin orientations 
(two out of three directions), which is either closed or ex- 
tends to infinity. All sites neighboring the loop share the 
common third spin direction. It is then possible to col- 
lectively rotate the spins on the loop freely around the 
third direction at zero energy cost. Such a new state is 
clearly non-coplanar, but still a ground state. 

Let us now discuss what carries over of this large clas- 
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FIG. 15: Classical GS's of the Heisenberg model on the 
kagome AFM which are selected by quantum or thermal fluc- 
tuations, (a) In zero-field, these are the 120° states. In finite 
fields ((b) and (c)), the three spins lie on the plane of the field, 
(b) For H < Hs/3, we have a one-parameter (0 < 120°) fam- 
ily of states with one of the spins (C) pointing antiparallel to 
the field, (c) For H > Hs/3, we have a two-parameter (^ and 
r;) family of states with two spins (A and B) being coUinear. 




(a)r(q=0) 



(b)M 



FIG. 16: Two of the 24 possible vertex three-colorings on the 
cuboctahedron which is projected on a plane here (periodic 
boundary conditions along the two arrows are implied), (a) 
One of the 6 colorings of the F (or q — 0) family, and (b) one 
of the 18 colorings of the Af family. 



sical degeneracy on the two present molecules. It is 
straightforward to enumerate all vertex-three colorings 
for both the cuboctahedron and the icosidodecahedron 
and the respective numbers are 24 and 60, or 4 and 10 
if one discards global recolorings. Two such states of the 
cuboctahedron are illustrated in Fig. 1161 while a typical 
one for the icosidodecahedron can be found in Fig. 2 of 
Ref.|43. We stress here that these 4 and 10 states are un- 
related by the global 0(3) symmetry, and therefore form 
genuinely different GS's and give rise to distinct "tower of 
states" at low-energies. The spatial symmetry properties 
of these states have been derived in Appendix [B] and can 
be summarized as follows. The 24 vertex three-colorings 
of the cuboctahedron form two invariant (under the op- 
erations of Ofi) families Pabc ^^"^ Pabc which consist 
of 6 and 18 states respectively, and they decompose into 
IR's of Oh as 



^ABC 

■pM 



Alg ® A2g © 2Eg , 
3(Alg©Eg©T2u) 



(19) 



As mentioned above, these are exactly the IR's that ap- 
pear (with open symbols) in the spectra of the lower two 
panels of Fig.[T31 On the other hand, the 60 vertex three- 
colorings of the icosidodecahedron form two invariant 
(under the operations of l/j) families which are equiva- 
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TABLE V: Heisenberg point: Decomposition of semiclassical 
coplanar states into IR's of G = SU(2) x R (where R = 0^ 
or \h) up to 5* = 6. N denotes the number of IR's for each S 
sector of SU(2), and is equal to (25' + 1) times the number of 
coplanar states in each family divided by six (i.e., 4 and 10 for 
the cuboctahedron and the icosidodecahedron respectively). 
For the derivation see Appendix [B] 



s 


Cuboc. 

Pabc(1 X 6) 


N 


Cuboc. 

Pabc(3x6) 


N 


Icosi. 

Rabc(10 X 6) N 





Alg 


1 


Alg,Eg 


3 


Ag,Au,Fg,Fu 10 


1 


A2g,Eg 


3 


Alg,Eg, 
2T2u 


9 


Ag,Au,Fg,Fu, 
2(Hg,Hu) 30 
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Alg,2Eg 
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3(Alg,Eg), 
2T2u 


15 


Ag,Au,Fg,Fu, 
4(Hg,Hu) 50 
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Alg,2A2g,2Eg 
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3(Alg,Eg), 
4T2u 


21 


3(Ag,Au,Fg,Fu), 
4(Hg,Hu) 70 
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2Alg,A2g,3Eg 
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5(Alg,Eg), 
4T2u 


27 


3(Ag,Au,Fg,Fu), 
6(Hg,Hu) 90 


5 


Alg,2A2g,4Eg 


11 


5(Alg,Eg), 
6T2u 


33 


3(Ag,Au,Fg,Fu), 
8(Hg,Hu) 110 


6 


3Alg,2A2g,4Eg 


13 


7(Alg,Eg), 
6T2u 


39 


5(Ag,Au,Fg,Fu), 
8(Hg,Hu) 130 



lent to each other (cf. Appendix [B|) . Here, we shall treat 
these families collectively as a single one called Rabc- 
This decomposes into IR's of 1^ a: 



,60 



Rabc = 2(Ag ® Au e Fg ® Fu e 2Hg 2Hu) . (20) 

As to the non-coplanar GS's (in zero- field), it turns out 
that the icosidodecahedron has none, since all the alter- 
nating loops described above have maximal length (20), 
and therefore the rotation of the spins on the loop just 
changes the global spin plane, thus preserving the copla- 
narity. The cuboctahedron however has non-coplanar 
GS's, since the loops can have length shorter than eight, 
in agreement with previous studies^i^. 

When switching on quantum fluctuations on the 
kagome lattice it is known that the spin waves at har- 
monic order select the coplanar GS's over the non- 
coplanar ones (order-by-disorder effect), due to the larger 
number of soft modes of the former—. A complete lifting 
of the remaining (spatial) degeneracy is taking place at 
the level of anharmonic spin waves, whereby the single 
-\/3 X \/3 magnetically ordered state is selected^i^. On 
the other hand for the extreme quantum case of s = 1/2 
a number of numerical works clearly show the absence 
of any magnetic orderi^ii. So based on these conflicting 
results it is difficult to predict to which regime the inter- 
mediate values of spin will belong. 

We now give a symmetry analysis for the cubocta- 
hedron spectra at low magnetizations which suggests 
strongly that the low-lying excitations can be described 
in semiclassical terms at the harmonic spin-wave level. 
To this end, we should first emphasize that the non- 



coplanar GS's do not carry the above spatial symmetry 
because not all triangles share the same spin plane in 
these configurations and thus spatial operations gener- 
ally cannot relate different triangles. This means that 
the striking agreement between the spatial symmetry of 
the above low-energy spectra and Eq. (|19p is not just acci- 
dental. Of course there can be other states whose spatial 
decomposition can in principle contain some or all of the 
IR's of Eq. (fT9| (in fact such states will be examined be- 
low for finite magnetizations). More stringent evidence 
comes by comparing the full spatial-|-spin symmetry pat- 
tern of the exact spectra to that of the 120° states. The 
latter has been derived in Appendix [B] and the results 
are provided in Table IVl up to 5 = 6. A closer inspection 
of the s = 3/2 and 5/2 spectra (cf. lower two panels of 
Fig. fn)) shows a remarkable agreement: All lowest-energy 
levels (shown with open (black) symbols) which are below 
the levels shown with filled (red) symbols can be identi- 
fied in TableFVlwith the right combinations of spatial and 
spin representations and multiplicity. It is important to 
note here that, although these towers are severely split by 
quantum fiuctuations — in fact some IR's of TableFvKe.g. 
one Alg level at the S = sector) can be found slightly 
higher than the lowest (red) filled symbols — almost the 
entire set of levels contained in Table |V] are found below 
the filled symbols with no extra level appearing. Thus 
we believe the combined spatial-|-spin symmetry pattern 
of the spectra for small magnetizations is a characteristic 
fingerprint of the 120° semiclassical states. 

Finite-field case. — Here the ground state manifold is 
generally larger since the classical constraint Sa = H/2 J 
allows for non-coplanar configurations already at the 
level of a single triangle. It is known^^'^^'^- however 
that the coplanar states with the spins lying in the field- 
plane have the largest number of soft modes and thus 
must be selected by quantum or thermal fiuctuations. 
A subsequent selection which depends on the field takes 
place within the field-plane for the orientation of the spin 
triad^^. It turns out that the most relevant GS's in a fi- 
nite field are the ones depicted in Figs. [T5f b) and (c)i^ 
For H < HJ3 (cf. Fig. [I^b)) the relevant GS's on a 
triangle form a one-parameter family with one spin anti- 
parallel to the field. It is important to note that, apart 
from their difference in the directions of the three spins, 
these configurations are spatially indistinguishable from 
the 120° states of Fig. [TSf a): They are both 3-sublattice 
states and thus carry the same spatial multiplicity (i.e. 
24 for the cuboctahedron and 60 for the icosidodecahe- 
dron) and the same spatial symmetry (i.e., Eqs. (|19p and 
(PO)) ). Quite similarly, the "quasi-coUinear" configura- 
tions of Fig. fTST c). which can be selected for H > Hg/'i, 
have an "uud" spatial structure. Hence they share the 
same spatial multiplicity and symmetry properties with 
the M = Ms /a Ising GS's. Namely, there exist 9 and 36 
"quasi-coUinear" states for the cuboctahedron and the 
icosidodecahedron respectively, and their spatial symme- 
try is given already in Eqs. (O and (fT^ . We should note 
here that the set of IR's appearing in Eqs. ([5]) and ([13 
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form a subset of the ones appearing in Eqs. (|19p and ([20]) 
respectively.^^ 

According to the above, all types of semiclassical con- 
figurations of Fig. [15] have the same set of spatial repre- 
sentations, except the "quasi-collinear" states which do 
not contain the A2g representation: At large magnetiza- 
tions this level is pushed higher in energy and this seems 
to confirm that it does not belong to the relevant towers 
of states of the "quasi-collinear" classical states. More 
generally, the fact that the same set of spatial IR's ap- 
pear in the low-energy spectra at all magnetizations signi- 
fies that our previous semiclassical interpretation for the 
zero-field case carries over for finite fields as well. Again, 
one may ask for more stringent evidence by a compari- 
son of the combined spatial-f-spin symmetry properties 
of these finite-magnetization states. We have derived 
these combined symmetries (not shown here) following 
the same lines as in Appendix [Bl but it turns out that 
the small size of the cuboctahedron together with the se- 
vere splitting of the low-lying states do not allow for a 
straightforward and thus definite identification of the rel- 
evant towers of states as above. As we show below, this 
will be possible for the much simpler case of the s > 1/2 
XY model. 



s=5/2 Cuboctahedron, XY point 



B. 



s > 1/2 XY model 



Here we study the XY model (i.e. the Jz = limit of 
Eqs. ([1])-(I1])) on the cuboctahedron. The reason of doing 
this is twofold. First because, in contrast to the Heisen- 
berg point, the XY point exemplifies very evidently the 
core idea of the simultaneous appearance of several tow- 
ers of states due to the spatial degeneracy of the classical 
GS's. And second, to provide an additional interpreta- 
tion of the Heisenberg spectra, since these can be thought 
of as being adiabatically connected to the XY spectra but 
split by the quantum fluctuations introduced by Jz- 

We first consider the classical case with a magnetic 
field perpendicular to the xy-plane. As above, the XY 
Hamiltonian can be rewritten in terms of the three spins 
Si of each triangle A and its total spin Sa as (in units of 
gfiB = 1): 



n. 



XY 

classical 



J. 



^E(si.-E' 



'i_L 



2H 



(21) 



where S^ 



^Ax~ 



S^ and similarly for s^^ . The lead- 
ing term of Eq. (|2ip is minimized by taking Saj^ = on 
each triangle. The remaining terms require that we max- 



imize both s^j^ and Si 



with the constraint s^j^ -|- sj^ 



s{s + 1) (the balance between the two components is con- 
trolled by the ratio 2H/Jxy). In zero-field this gives the 
3-sublattice states where all triangles share the same spin 
(xy) plane. The major difference with our previous anal- 
ysis of the Heisenberg (classical) GS's is that here a spin 
plane (i.e. the xy plane) is selected explicitly from the 




0.1 



0.08- 



0.06 - 



0.04 



0.02 



s=5/2 Cuboctahedron, XY point 



1 


1 


1 1 1 1 1 1 1 1 1 

OXl <| 


- 


OAlg 
n A2g 
< Eg 

>T2u 




IVI family W 
r family © 


-pOO©©^CKl<Q^<0<]<p<l<10<l<pO<MO<^ 



10 



15 

s. 



20 



25 



30 



FIG. 17: (Color online) Low-energy spectra of the XY model 
on the s = 5/2 cuboctahedron. The arrow in the upper panel 
indicates the set of Anderson towers of states (or rotational 
bands). These are shown in finer energy resolution in the 
lower panel in order to reveal their symmetry structure. The 
latter is in full agreement with Tabl dVll Note the level cross- 
ing between the two semiclassical families occurring slightly 
below the M = Ms/3 field. 



beginning (i.e., for zero field) and this gives a finite en- 
ergy cost to non-coplanar configurations. A finite field 
gives rise to a tilt of the three spins out of the xy plane 
giving rise to the so-called "umbrella" states. It is clear 
that both in zero and in finite field, the classical XY GS's 
have the same spatial multiplicity and spatial symmetry 
properties (but not spin symmetry properties, see below) 
with the set of coplanar 3-sublattice GS's of the Heisen- 
berg point. 

Let us now consider the quantum-mechanical XY 
model. For our demonstration purposes, it suffices to 
consider the s = 5/2 case only (the s = 3/2 case is very 
similari^) . The low-energy spectrum is shown in the up- 
per panel of Fig. [T71 A number of spectral features are 
revealed. First, in contrast to the Heisenberg case studied 
above, the lowest-energy portion of the spectrum (indi- 
cated by the arrow) is well isolated from higher excita- 
tions and it comprises four distinct towers of states. This 
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TABLE VI: XY point: Decomposition of semiclassical copla- 
nar states into IR's of G = Ccxju x R (where R = 0^ or U). N is 
the number of IR's for each sector of Ccx)ti and is equal to the 
dimensionality of the corresponding IR of Good fcf.Table lVIIip 
times the number of coplanar states in each family divided by 
six (i.e., 1 for Pabc, 3 for Pabc. and 10 for Rabc)- The 
details of the derivation are given in Appendix iBl 
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Cuboc. 

PiBc(l X 6) N 


Cuboc. 

PfBc(3 X 6) N 


Icosi. 
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Alg 1 
A2g 1 
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2Hg,2Hu 20 
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Alg,A2g 2 
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multiplicity is a fingerprint of the spatial degeneracy of 
the classical 3-sublattice states mentioned above. This 
can be further substantiated by examining more closely 
the symmetry structure of the excitations intervening in 
these towers. To this end we zoom in on these towers 
in the lower panel of Fig. 1171 A simple inspection of the 
spatial IR's that appear in this panel reveals that they 
are exactly the ones given by Eq. ([19]). Much stronger 
evidence comes by examining the full spatial+spin sym- 
metry structure of the lowest towers. Indeed, a closer 
comparison to Table IVTl (derived in Appendix [BJ demon- 
strates that there is a remarkable one-to-one correspon- 
dence of each of the lowest towers with the classical fami- 
lies which holds in almost the entire magnetization range. 
We should emphasize here that each semiclassical state 
shows a different and quite non-trivial symmetry pattern 
which is in some sense a very characteristic fingerprint 
of the state. For instance, the full content of 24 states 
of Eq. p^ is recovered every three Sz sectors in the 
lowest towers with the specific pattern of combined spa- 
tial and spin IR's given in Table IVII This remarkable 
agreement between exact ED spectra and our symme- 
try derivation is indeed a strong evidence that the lowest 
towers of states can be thought of as renormalized semi- 
classical 3-sublattice states. 

Some additional remarks are in order here regarding 
the energies of the two families of towers as revealed in 
the lower panel of Fig. [T71 We should first note that all 
towers are expected to become degenerate in the classical 
s ^ 1 limit. Based on spectra with s = 3/2 (not shown 
here), we find that the energy splittings within a tower 
of each of the two families P^^^ and Pabc diminishes 
quickly on increasing s as expected, while at the same 
time, the energy splitting between them remains sizable. 
On the other hand, as a function of Sz (or field), there 
is an interesting level crossing between the two families 
somewhat below the M = Ms/3 magnetization plateau, 
with the M family being more favorable below this point. 
An understanding of the above level-crossing could arise 
by employing for instance a semiclassical expansion for 
the XY model in a field, in a similar fashion with what is 



done for the Heisenberg model^i^, but such an analysis 
is clearly beyond the scope of this paper. 



V. SUMMARY 

We have presented an extended study of the low-energy 
physics of two existing magnetic molecule realizations 
of the kagome AFM on the sphere, the cuboctahedron 
and the icosidodecahedron. Our ED results revealed 
a number of generic spectral features which stem from 
the corner-sharing topology of these clusters. Indeed, a 
simple comparison to a finite-size s = 1/2 unfrustrated 
magnet demonstrated that frustrated clusters manifest 
a "bulk" of very dense low-energy excitations. We fo- 
cused on two major aspects which are of general inter- 
est but were particularly oriented toward the s = 5/2 
Mo72Fe3o cluster: (i) the low-energy excitations above 
the M = Ms/3 plateau and (ii) the low-lying spectra of 
the Heisenberg model for s > 1/2. 

For the M = Ms/3 plateau, we first demonstrated that 
the s — 1/2 icosidodecahedron shows 36 low-lying ex- 
citations which are adiabatically connected to collinear 
"uud" Ising (GS's), at the same time being well isolated 
from higher levels by a relatively large energy gap. We 
then argued, based on a complementary physical picture 
which emerged from the derivation of an effective quan- 
tum dimer model, that this s — 1/2 feature must be spe- 
cial to the topology of the icosidodecahedron and that 
it must survive for s — 5/2 as well. We also predicted 
that the corresponding 36 low-lying plateau states of the 
s = 5/2 icosidodecahedron consist of two "uud" families 
(of 30 and 6 states respectively) which are separated by a 
small diagonal energy splitting. This result can be con- 
firmed by low-temperature specific heat measurements 
at the M = Ms/3 regime {H ~ 5.9 Tesla) of Mo72Fe3o 
and/or by an assessment of the associated missing en- 
tropy. 

In the second part, we showed exact diagonalization 
spectra for the s > 1/2 Heisenberg cuboctahedron which 
demonstrated that the dense low-lying excitation fea- 
tures of the s = 1/2 case are present for s > 1/2 as 
well, albeit with a striking spatial-|-spin symmetry pat- 
tern. These spectra provide a semiclassical interpreta- 
tion of the broad inelastic neutron scattering response 
reported for Mo72Fe3o. The main ingredient of this in- 
terpretation is the simultaneous presence of several low- 
energy towers of states or rotational bands at low energies 
which originate from the large spatial multiplicity of the 
classical Heisenberg ground states, and this is known to 
be a generic feature of highly frustrated clusters. This 
semiclassical interpretation was further corroborated by 
an independent group theoretical analysis which demon- 
strated that the striking symmetry pattern of the low- 
lying excitations is indeed a characteristic fingerprint of 
the classical coplanar ground states. The core idea of the 
simultaneous presence of several rotational bands at low 
energies due to the discrete classical degeneracy was fi- 
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nally exemplified very evidently by a study of the s > 1/2 
XY model. 



VI. ACKNOWLEDGMENTS 

We would like to acknowledge fruitful discussions with 
C. L. Henley, M. Luban, and K. P. Schmidt and ear- 
lier collaboration with J.-B. Fouet and S. Dommange 
on related subjects. The point group symmetry data 
used in our calculations have been obtained using the 
Bethe package^. We would like to thank K. Rykhlin- 
skaia for her assistance with this package. This work was 
supported by the Swiss National Fund and by MaNEP. 
The computations have been enabled by the allocation 
of computational resources on the machines of the CSCS 
in Manno. 



APPENDIX A: DEGENERATE PERTURBATION 
THEORY AROUND THE ISING POINT 

Here we describe some very general considerations 
which greatly facilitate the classification of processes ap- 
pearing in degenerate perturbation theory around the 
Ising point. We consider corner sharing triangles struc- 
tures with general spin s, and focus on the M = Ms/3 
plateau. The ground-state Ising manifold consists of con- 
figurations which maximize the number of extremum lo- 
cal moments m^". At the plateau M = Ms/3 phase, 
this consists of all configurations with two spins having 
m — s and one with ttt, = — s in each triangle. All pro- 
cesses triggered by Tixy must preserve this constraint. 
Let us denote by Eq the zero-th order energy and by 
'Pq,Qo = 1 — 'Po the projections onto and out of the Ising 
manifold. We also designate by TZ the resolvent operator 



n ^ (Eo - QoHoQo)- 



(Al) 



The n-th order term of the effective Hamiltonian in the 
Rayleigh-Schrodinger formulation^ reads 

H^ff^ = Vonxy{nn.j:yy'-^Vo+ remaining terms , (A2) 

where each of the "remaining terms" can be thought of 
as a product combination of lower order terms n^^f (with 
k < n). This separation is useful when e.g. all processes 
below some order n are constant since then, to order n, 
it suffices to keep only the leading term of Eq. (|A2p . The 
presence of Pq's in Eq. (jA2|) enforces all terms to flip 
spins in such a way as to respect the "uud" constraint in 
each triangular unit. 

The derivation of Tioff at any given order is greatly fa- 
cilitated by noting that only "linked" processes should 
be taken into account. These are interactions that are 
"connected" in the following sense. Substituting Hxy in 
Eq. (jA2p gives different types of terms, each one carrying 
a string of a given number of bond operators s~ sj' . We 
differentiate between "linked" and "unlinked" interaction 







(c) 



FIG. 18: As above, filled circles denote spins with m = —s, re- 
maining vertices have m — s. (a) All configurations of a single 
triangle have the same energy to all orders, (b) The configu- 
rational energy of two adjacent triangles depends only on the 
shared spin, since all 4 states with the shared spin pointing up 
have the same energy, (c) The contraction method is based 
on the observation that E{R; si, S2, sg) — E{R; si). 



terms depending on whether the set of all vertices appear- 
ing in the corresponding string forms a connected (open 
or closed) path in the lattice or not. The latter contain 
nonlocal interactions between disconnected parts of the 
lattice and therefore must be omitted^^ . Only connected 
paths should therefore be considered. Further simplifica- 
tions arise from the "uud" constraint as described below 
for diagonal and off-diagonal processes. 

Diagonal processes. — As explained in considerable de- 
tail in Ref. KJ in the context of the pyrochlore lattice, all 
diagonal processes up to a given order give an overall con- 
stant energy shift, with the leading non constant terms 
arising from processes along closed loop configurations. 
As we explain below, similar results apply to the present 
case of corner sharing triangles as well. The proof can be 
demonstrated in a compact way by using the contraction 
method of Bergman et alJ^ 

Consider first all-order diagonal processes confined to a 
given triangle. Clearly, the only physically distinct con- 
figuration on this triangle is the "uud" one, since the 
associated diagonal energy does not depend on which of 
the three vertices the down spin resides (cf. Fig. [TST a)). 
By sampling the energies of all triangles we end up with 
an energy that is global, i.e., the same for all Ising states. 
Hence all-order processes confined in a single triangle 
give an overall constant energy shift and can thus be 
neglected. We now consider all-order processes confined 
to two adjacent triangles. Here, there are two physically 
distinct classes of configurations (cf. Fig.fTSTb)). depend- 
ing on whether the shared spin points down (first state in 
Fig. [T8jb)) or up (remaining four states). Since all four 
states with the shared spin pointing up are physically 
equivalent, the energy is only a function of the shared 
spin variable. By sampling again over the lattice we ob- 
tain a global energy shift. Remarkably, these arguments 
can be generalized to much broader cases by the contrac- 
tion method exemplified in Fig. llSf c): Since permuting 
the spin variables 2 and 3 results in a topologically equiv- 
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FIG. 19: A fragment of the corner- sharing triangle lattice 
with no closed loops. This is a Bethe lattice made of triangles, 
known as Husimi cactus. 



(a) 



(b) 
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FIG. 20: The thick solid line denotes a path corresponding 
to a given term in the perturbation series. Only the alternat- 
ing spin up-down configuration shown in (b) (which must be 
closed) is amenable to an ofi^-diagonal kinetic process. 



alent state, knowing the value of si is enough, i.e., 

E{R;si,S2,S3)=E{R;si) , (A3) 

where R designates all remaining spin variables. This 
"contraction" can be continued with the next available 
triangle inside the shaded area R, and so on. If this pro- 
cess can be continued until we are left with a function of 
a single vertex, then sampling this function over all spins 
we obtain again a global energy shift. A non-constant 
energy contribution may arise only when the contraction 
process cannot be continued until the last spin variable. 
This happens whenever the shaded area R of Fig. \TE[ c) 
contains one or more closed loops, since none of the tri- 
angles making up a loop is "contractible" . This leaves 
us with the following quite general statement: All-order 
processes confined to a fragment of the lattice with no 
closed loops give an overall constant energy shift. The 
most general form of such fragments is depicted in Fig.fTOl 
and is recognized to be a Bethe lattice made of triangles, 
known as Husimi cactus. 

Hence, the lowest order diagonal processes come from 
closed loops in the lattice. One should also remark that, 
in contrast to off-diagonal processes (see below), the or- 
der at which diagonal processes first appear is indepen- 
dent of the spin s. 

Off-diagonal processes. — Figure HOI shows schemati- 
cally a given lattice path (solid thick line) with two par- 
ticular configuration choices. All spins have m = s except 
the ones indicated by filled circles with m = —s. By defi- 
nition, only spins residing on this path may be flipped. In 



TABLE VII: Degenerate perturbation theory calculations for 
the off-diagonal kinetic amplitude t on alternating up-down 
configurations around loops with L — 4, 6 and 10 sites, and 



for 


various 


intrinsic spins 


s. The order in 


Hxy and the total 


number of 


contributing processes are also 


given. 


L 


2s 


Order (is) 


t [J.a^^] 


# of processes 


4 


1 


2 


-1.0 


4 




2 


4 


-1.0 


36 




3 


6 


-0.5625 


400 




4 


8 


-0.25 


4 900 




5 


10 


-0.09765625 


63 504 




6 


12 


-0.03515625 


853 776 


6 


1 


3 


1.5 


12 




2 


6 


-0.88402469 


900 




3 


9 


0.25093125 


94 080 




4 


12 


-0.05637473 


11 988 900 




5 


15 


0.01106939 


1 704 214 512 




6 


18 


-0.00199964 


260 453 217 024 


8 


1 


4 


-2.5 


48 




2 


8 


-0.93709194 


45 360 




3 


12 


-0.12770306 


60 614 400 




4 


16 


-0.01464521 


114 144 030 000 


10 


1 


5 


+4.375 


240 




2 


10 


-1.0924858 


3 855 600 



off-diagonal processes at least one spin, say Sq, is flipped 
from m — ±s to to = =ps. But, since the final state must 
preserve the "uud" constraint in each triangle, flipping 
So must be accommodated by a similar flip of the adja- 
cent spins lying on the path, namely s_i and Si. Simi- 
larly, the spins at vertices 2 and 3 must also be flipped. 
If, as drawn in Fig. [20y a). the initial conflguration has 
TO.3 — TO,4 = — s, then flipping 3 already violates the 
"uud" constraint. It is also clear that if the path is open 
at one end (or at both ends) the ending spin (or spins) 
will be flipped at the flnal state, thus violating the "uud" 
constraint on the adjacent triangle(s). Thus only closed, 
connected paths with alternating up-down spins, such as 
in Fig. [SOT b) are amenable to an off-diagonal process. 
The simplest possible paths are loops of even number of 
spins L. Since flipping each spin requires 2s operations, 
off-diagonal processes on a simple loop appear in order 
a^^. An example was shown in Fig. [9] for the icosidodec- 
ahedron. In Table IVIII we provide calculated numerical 
values for the amplitude of such virtual processes around 
loops with L = 4, 6, 8 and 10 sites and various values of 
s. Starting from a given loop, one may also build pro- 
cesses (and paths) of higher order by invoking adjacent 
triangular units. An example for the icosidodecahedron 
was the fifth order processes mentioned in Sec. IIIIB21 
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APPENDIX B: SYMMETRY PROPERTIES OF 
THE 3-SUBLATTICE COPLANAR STATES 

Here we give the details of the derivation of the symme- 
try properties of the semiclassical 3-sublattice coplanar 
states of Fig. fTST a). These are relevant for the Heisen- 
berg model at iJ = as well as the XY model (dis- 
cussed in Sec. lIVSp for both zero and finite fields. As 
explained above, the spatial symmetry properties of the 
states of Fig. [TST b) are the same as that of (a), and simi- 
larly the spatial symmetry of the "quasi-coUinear" states 
of Fig. [Wc) is identical to that of the "uud" states given 
in Eqs. (O and ([T^ . The derivation of the combined 
spatial+spin properties of Fig. [TSf b) and (c) are not of 
our interest here but can be found easily following the 
same steps as below. 

We are using the following notation and conventions. 
The groups of real space and spin space operations are 
denoted respectively by R and L. In particular, R = 0^ 
for the cuboctahedron and l/j for the icosidodecahedron. 
The full group R x L is designated by G. A stabilizer 
He of a classical state |c) consists of elements h which 
preserve |c) (i.e. h\c) — \c)) and this will be either a 
subgroup of R or a subgroup of G depending on whether 
we are examining only the spatial or the full symmetry 
properties. The elements of R, L and G are labeled by 
r, I, g = r ■ I respectively, while their IR's are denoted 
as T>P{r), D^(Z), and B^g) = D''(r) ® D^(0. Similarly, 
their characters are denoted by x''(r), x^(0; ^^'^ x'ig = 
r ■ I) = x^if) • X^iO- Let us first discuss the spatial 
symmetry and then the combined spatial+spin symmetry 
structure of the 3-sublattice states. 



1. Spatial symmetry of 3-sublattice states 

As we discussed previously, each vertex three-coloring 
|c) is in one-to-one correspondence with the states of 
Fig- fW aV Starting from a given \c) and applying all 
elements r of R we generate an invariant vector space or 
orbit O. The decomposition of O into IR's D^''' of R is 
given by the well known formula^ 



0(r) = ^mpD''(r), 
p 

"^P = ME^''W*Tr[0(r 



(Bl) 

(B2) 



reR 



The matrix element Occ'{r) — {c\r\c') is equal to one if r 
belongs to the stabilizer He of \c) and vanishes otherwise. 
Thus we may rewrite Eq. (jB2p as 



(B3) 



^m,S x'C). 



heUaCR 



where we made use of |R| = |0| • |Hc|. This relation 
follows from the coset decomposition of R with respect 
to He (each coset is in one-to-one correspondence with 



the states |c') of the orbit O). Employing (jBSp to all 
different orbits we obtain the symmetry properties of the 
coplanar states. 

Let us see now what happens for the cuboctahe- 
dron and the icosidodecahedron separately. As discussed 
above, the cuboctahedron has a total number of 24 vertex 
three-colorings \c) = |ABC). Under Oh they form four in- 
variant orbits of six colorings each. The first orbit, called 
■f*ABC consists of the six global permutations of the trans- 
lationally invariant coloring depicted in Fig. I16f a). The 
18 colorings that belong to the remaining three orbits 
result from the first orbit by interchanging colors along 
loops with two alternating colors. One such configuration 
is shown in Fig. [TBFb). Although these three orbits are 
equivalent it is useful for the following discussion of the 
full spatial-f spin properties to treat them collectively as 
a single one which we term Pabc'''^ "^^^ decomposition 



of the above orbits into IR's of 0^ was given in Eq. (|19p . 
On the other hand the icosidodecahedron has 60 color- 
ing states. Under \h they form two orbits of 30 colorings 
each. The first orbit consists of the 3 cyclic permuta- 
tions (ABC,CAB,BCA) of 10 ABC states, while the sec- 
ond orbit consists of the remaining 3 permutations (or 
reflections) (BAC,CBA,ACB) of the same 10 states. Al- 
though these two orbits are equivalent it is useful for the 
discussion of the full spatial-|-spin properties (see below) 
to treat them as a single one which we denote by Rabc- 
Its decomposition into IR's of l/j was given in Eq. ([20|l . 



2. Spatial+spin symmetry of 3-sublattice states 

We shall now go one step further and derive the 
combined spatial+spin properties of the above states. 
Namely, a decomposition similar to that of Eqs. p^ and 
([20|) but now in terms of IR's D'*' of the full symmetry 
group G = R X L of the Hamiltonian. The method has 
been employed previously in the seminal works of Bernu 
et flZ;. S .-i i. i - i and Lecheminant et al^'^. 

The recipe is quite analogous to the one we employed 
above. Here however, by applying the elements of the 
full group G on a given classical state we generate a con- 
tinuous orbit O. All states contained in this orbit are 
coplanar but their order parameter has all possible orien- 
tations in spin space (O is a continuous "order parameter 
space" ) . Otherwise, the equation giving the numbers m^ 
(i.e. how many times is D"*" appearing in the decomposi- 
tion of O into IR's of G) is fully analogous to Eq. (jB3p 
and reads 



m^ 



1 

Wr. 



E ^'^^)^ 



(B4) 



/leHeCG 



where now x'^i^ — r ■ I) = %''(?') • x'*'(0- For the identifi- 
cation of He it is expedient to split the operations of the 
space group R into the 3! sets Sabc, Scab, • ■ • , ^cba defined 
as follows. The first consists of elements which map \c) 
to itself. On the other hand, the elements of Seah rnap 
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TABLE VIIL Character table of Coo« (see e.g. Ref.HJ). Here, 
E denotes the identity, R^ the set of U(l) rotations, and cr„ 
the set of aU vertical mirror planes. The IR's Coou are char- 
acterized by \Sz\, and the "parity" under (t„ for S^ = 0. 



15^1, a„ 


E 




R4) 


cr„ 


0,+ 


1 




1 


1 


0,- 


1 




1 


-1 


n> 1 


2 


2 


cos 71 (^ 






TABLE IX: The numerical values of the integers Nabc, Ncab, 
etc. defined in the text. Remaining IR's or blank entries 
correspond to vanishing values. We also give the order of the 
corresponding stabilizers He C G. 





Cuboc. r 


Cuboc. M 


Icosi. 




(Alg,A2g,Eg) 


(Alg,Eg,T2u) 


(Ag,Au,Fg,Fu,Hg,Hu) 


|Hc| 


48 


16 


12 


Nabc 


(8, 8,16) 


(8, 8, 8) 


(4, 4, 4, 4, 8, 8) 


Ncab 


(8, 8,-8) 




(4, 4, 4, 4,-4,-4) 


Nbca 


(8, 8,-8) 




(4, 4, 4, 4,-4,-4) 


Nbac 


(8,-8, 0) 


(8, 8,-8) 




Nacb 


(8,-8, 0) 






Ncba 


(8,-8, 0) 







\c) to its globally permuted (ABC) 1-^ (CAB) version, 
and similarly for the remaining sets. We also define the 
following set of integer numbers 



N' 



be 



heSab 



X'{h), <, 



xP{h), etc . (B5) 



heSe 



These numbers depend on the transformation properties 
of \c) under the spatial group R alone. The non- vanishing 
ones are given in Table IIXI 

We should note here that one can make a choice of G 
depending on the amount of information we seek. For 
instance we may choose according to the symmetries we 
implement in our exact diagonalizations. We may even 
take L as the idenity, i.e. G = R. In the latter case we re- 
cover Eq. (|B3p . The type and number of invariant vector 
spaces in each case will be different and the symmetry 
decomposition must be applied to each one separately, 
but the corresponding results will be consistent with each 
other. 

Let us now apply the above to the zero-field Heisenberg 
model and the XY model. 



a, SU(2) point 

Here we take L = SU(2) and A is the total spin S 
which is integer here. A single |c) generates the full set 
of coplanar states for each family of the clusters. The 
elements of Sbac, ^cba, and Sacb can be combined with 



TV rotations of SU(2) and bring \c) back to itself. Thus 

again |Hc| = \Sabc\ + h |Sc6a|- Using Eq. JBU with 

7 = (p, S) then 



rrip^s 



KbcX^'Ho) + KabX^^Hf) + <.x(^Hf ) 

+ «c + Kcb + ^L) X^'\^)] /|Hc| . (B6) 



where^ X^^H't') = ^'"in(f>/2 • Replacing the values of 
Table IIXI for each family separately we obtain the sym- 
metry structures given in Table [V] 



b. XY point 

Here, we take L = Coo-u which includes U(l) rotations 
and the continuous set of vertical (i.e. containing the z- 
axis) mirror planes. The IR's of Good can be generally 
labeled'^^ by a non-negative integer n or \Sz\, with an 
additional label cr„ = ±1 for n = which stands for the 
parity under the mirror operation. Hence A = (|S'2|,cr^). 
All IR's for n = \Sz\ ^0 are two-dimensional and consist 
of pairs of Sz and —Sz basis vectors. The characters of 
Good are given in Table IVini 

It suffices to select a single coloring state \c) since this 
generates all coloring states for both clusters. The stabi- 
lizer He can be found as follows. Each one of the sets Sabc, 
Scab, etc. discussed above can be combined with one of 
the elements of Goo« to give |c) again. For instance, an el- 
ement of Scab can be combined with a U(l) spin rotation 
of 2tt/3. On the other hand, an element of Sbac can be 
combined with a mirror plane containing the C-axis (i.e. 
that of the spins "colored" as C) . The set of all such com- 
bined operations span He, i.e. |He| = \Sabc\ + • . - + \Sbac\- 
Using Eq. (|B4)) with 7 = (p. A) then 



mp^x 



[Kbcx\E) + n:^,x\^) + nLx\y' 



<N, 



P 

bac 



Kcb 



<.a)x'K)]/|H. 



(B7) 



Replacing the values of Table IIXI and the characters x^ 
from Table lylTTl for each family separately we obtain the 
symmetry structures given in Table |VT1 According to this 
table, the symmetry pattern repeats itself every three 
Sz sectors. In particular, the full set of spatial IR's of 
Eqs. (|19p and (PO)) is contained in any triad of subsequent 
Sz sectors. As can be seen from the above relations, 
this periodic pattern stems from the 120° 3-sublattice 
symmetry structure of the coplanar states (the character 
X^{4') ^ cos {n(f>) (with = 0, 27r/3,47r/3) has a period 
of n = 3). 

In connection to a remark above, it is clear that we 
could have chosen here L — U(l) instead of Goou- The 
corresponding stabilizer would then obviously be differ- 
ent from the above since none of the combinations of Sbac, 
Sacb and Scba with U(l) rotations can bring \c) to itself. 
Nevertheless, the results from the two different choices of 
L are consistent with each other. 
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We should finally emphasize a non-trivial feature which 
appears in both Tables IVTl and IVl and holds for each fam- 
ily separately. Namely that the total number Tn\ of states 
(counting the degeneracy dp of spatial IR's) for a given 
IR A of the spin group L is equal to the dimensional- 
ity d\ of this IR (i.e. (25* +1) for the Heisenberg case, 



see Table Em] for the XY case) times the ratio |R|/|Hc|. 
This feature can be proven by group theory alone using 
m\ = J2n^p "^p.A, Eq- (|B4|1 and the so-called character 
orthogonality relational. 
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